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1.  Introduction 

The  coronal  temperature,  often  exceeding  one  million  degrees,  is  a  direct  manifestation  of 
coronal  heating  processes.  While  this  remarkable  phenomenon  remains  to  a  great  extent 
unaccounted  for  by  the  physical  processes  that  we  are  currently  aware  of,  it  has  fueled 
endless  debates  and  countless  research  projects. 

We  have  attempted  to  explore  some  aspects  of  this  problem  by  pursuing  two  complementary 
directions:  (1)  the  improved  knowledge  of  the  temperature  structure  of  the  solar  corona, 
and  (2)  the  study  of  the  morphology  of  the  magnetic  field  that  defines  it.  Our  approach 
has  involved  data  acquisition  of  multiwavelength  coordinated  observations,  image  processing 
and  data  analysis. 

This  final  report  summarizes  our  main  research  projects  and  findings  which  have  been  accom¬ 
plished  with  support  from  the  Air  Force  for  Scientific  Research  from  July  1991  to  September 
1994.  A  list  of  publications  resulting  from  these  studies  is  given  in  Section  5. 

2.  Summary  of  Results 

Our  approach  has  focused  on  exploring  the  physical  characteristics  of  the  coronal  heating 
mechanisms,  as  manifested  within  different  large  scale  magnetic  structures,  such  as  coronal 
holes,  quiet  regions  and  active  regions.  These  explorations  were  made  using  different  data 
sets,  data  analysis  techniques  and  image  processing  tools.  We  present  first  the  characteristics 
of  the  data  used  (section  2.1),  we  then  describe  the  different  analyses  and  image  processing 
tools  developed  to  study  the  plasma  conditions,  such  as  temperature,  the  characteristics  of 
the  variable  emission,  and  the  morphology  of  the  open  and  closed  structures  in  the  solar 
corona  (section  2.2).  We  highlight  the  results  of  these  studies  in  section  2.3. 

2.1  Data  Characteristics 

2.1.1  Space-based  Observations 

One  of  the  outstanding  characteristics  of  the  EUV  Skylab  data  obtained  with  the  Harvard/S- 
055  Extreme  Ultraviolet  (EUV)  spectroheliometer,  was  that  the  spectroheliograms  were 
made  simultaneously,  through  the  same  instrumental  slit,  in  several  wavelengths  characteris¬ 
tic  of  chromospheric  to  coronal  emission.  To  date,  these  characteristics  have  been  unmatched 
by  any  subsequent  mission.  The  field  of  view  of  the  instrument  was  typically  5’  x  5’,  except 
during  special  observing  sequences  when  smaller  fields  of  view  were  scanned  to  allow  for  a 
higher  temporal  resolution  (see  Reeves  et  al.  1977  a,b).  The  field  of  view  was  rastered  every 
5.5  minutes  with  a  spatial  resolution  of  5”  x  5” .  Observations  were  made  with  different  spec¬ 
trometer  grating  settings.  In  the  standard  position,  Mg  X  625  A  (1.12  x  106K),  O  VI  1032  A 
(2.85  x  105K),  O  IV  554  A  (1.5  x  105K),  C  III  977  A  (7.4  x  104K),  C  II  1336  A  (3.16  x  104K), 
and  Lyman  a  1216  A  (2  x  104K)  were  recorded.  Alternative  grating  positions  yielded,  for 
example,  Si  XII  521  A  (1.8  x  106K)  and  Ne  VII  465  A  (5  x  105K),  together  with  other  lines 
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sensitive  to  chromospheric  emission.  (The  temperatures  listed  in  parentheses  correspond  to 
the  peak  temperature  of  formation  of  the  corresponding  lines). 

While  our  work  has  been  based  primarily  on  the  EUV  Skylab  data,  we  have  also  explored 
selected  observations  from  Yohkoh.  The  Yohkoh  spacecraft  was  launched  in  1991,  and  since, 
has  returned  x-ray  observations  of  the  Sun  with  an  unprecedented  temporal  coverage  (see 
Tsuneta  et  al  1991,  Acton  et  al  1992).  For  our  research  we  have  been  using  the  observations 
obtained  with  the  soft  x-ray  grazing  incidence  telescope  (SXT)  onboard  Yohkoh  which  forms 
whole  sun  images  in  the  wavelength  range  2-50  A  with  a  pixel  size  of  approximately  5”. 
Images  of  selected  regions  of  the  Sun  can  also  be  obtained  with  a  higher  spatial  resolution 
and  temporal  cadence. 


2.1.2  Ground-based  Observations 

Our  ground-based  observations  consist  primarily  of  data  from  the  National  Solar  Observatory 
at  Sacramento  Peak  (NSO/SP),  and  the  Very  Large  Array  (VLA).  Observations  of  the  Fe  X 
6374  A  and  Fe  XIV  5303  A  lines  were  made  with  the  Fisher-Smartt  Emission  Line  Coronal 
Photometer  (ELCP)  (Smartt,  1982).  This  instrument  photoelectrically  records  the  solar 
corona  when  fed  with  the  Evans  Solar  Facility  40  cm  Coronagraph.  Scans  are  recorded  around 
the  solar  limb  with  a  1.1’  aperture,  in  angular  steps  of  3°  in  latitude.  The  coronagraph  allows 
observations  to  start  at  a  heliocentric  distance  of  1.1  Rs  and  extend  reliably  out  to  1.2-1. 3 
Rs.  Photographic  images  of  the  corona  in  H  a,  Fe  X  and  Fe  XIV  were  also  obtained  with 
the  one  shot  coronagraph.  Image  processing  techniques  were  applied  to  the  subsequently 
digitized  images. 

Observations  were  also  obtained  at  several  cm  wavelengths  with  the  radio  telescopes  at  the 
Very  Large  Array  (VLA),  in  coordination  with  other  ground-  (e.g.  NSO/SP)  and  space- 
based  observations  (e.g.  SPARTAN,  Yohkoh ,  and  NIXT).  With  this  instrument,  the  field  of 
view  depends  on  the  wavelength  used,  and  the  spatial  resolution  depends  on  both  wavelength 
and  configuration  of  the  radio  antennas.  Our  observations  were  made  typically  at  6,  20  and 
90  cm. 


2.2  Data  Analysis  Tools 
2.2.1  Line  Ratio  Technique 

The  inference  of  the  coronal  temperature  from  the  ratio  of  the  intensities  of  emission  lines 
depends  on  the  assumption  that  the  density  and  temperature  are  constant  along  a  given  line 
of  sight  from  which  the  emission  from  different  spectral  lines  is  observed.  In  the  absence 
of  direct  measurements  of  the  coronal  temperature,  the  line  ratio  technique  is  the  most 
commonly  used.  To  apply  this  technique,  one  compares  the  ratios  of  line  intensities  deduced 
from  observations  with  those  derived  theoretically.  The  intensity  of  an  emission  line  at  a 
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given  wavelength  A  can  be  written  as 

=  [  P\NeNHds ,  (1) 

in  er#s  cm"2  where  the  emissivity,  PA,  is  given  by  (see,  e.g.  Raymond  1988,  and 

references  therein) 

ft  =  8.63  10-  ^  ^  hv  .  (2) 

In  these  expressions,  the  integral  is  along  the  line  of  sight,  L.  Ne  and  NH  are  the  electron  and 
hydrogen  (proton)  densities  respectively.  The  elemental  abundance  Ael  is  given  with  respect 
to  hydrogen,  ^  is  the  ionization  fraction,  Q,  is  the  collision  strength,  u  is  the  degeneracy 
of  the  ground  state,  and  E  is  the  threshold  energy  usually  taken  to  be  equal  to  the  photon 
energy  hu.  It  is  clear  from  this  expression  that  the  ratio  of  the  emissivity  (or  equivalently, 
intensity)  of  two  lines  from  different  elements  then  reduces  to  a  functional  dependence  on 
temperature,  as  well  as  on  the  ratio  of  the  abundances  of  the  ions  involved.  Since  the 
elemental  abundances  are  poorly  known,  no  unique  temperature  can  be  determined  from 
only  one  line  ratio.  However,  with  the  availability  of  three  or  more  lines,  and  hence  three 
or  more  ratios,  one  can  restrict  the  temperature  range  thus  inferred  (see  Habbal,  Esser  and 
Arndt  1993).  The  uncertainties  are  significantly  alleviated  when  different  lines  from  the 
same  element  can  be  measured. 

2.2.2  Image  Difference  Technique 

To  characterize  the  properties  of  the  intermittent  component  of  the  coronal  heating  processes, 
independently  of  the  magnetic  structures  studied,  i.e.  whether  in  coronal  holes,  quiet  or 
active  regions,  we  expanded  on  a  method  developed  by  Habbal  and  Grace  (1991).  The 
approach  is  as  follows:  the  change  in  the  intensity  in  any  of  the  pixels  making  up  a  given 
image  raster  is  considered  significantly  variable  if  the  absolute  value  of  the  difference  in 
intensity  in  a  given  pixel,  at  two  consecutive  times,  is  greater  than  2  <j,  i.e.  if 

51  =  |4  -4|  >  2a  (3) 

where  a  =  vT,  /  =  PP—2.  is  the  average  of  the  intensities  Itl  and  It2  within  each  pixel 
of  the  images,  measured  at  times  t\  and  Tz.  y  then  represents  the  relative  change  of  the 
intensity  between  two  consecutive  time  frames  for  which  the  variability  is  significant.  We 
also  introduce  the  measure,  y,  or  spatial  density  of  variability,  defined  as  the  ratio  of  the 
number  of  variable  pixels  satisfying  condition  (1)  to  the  total  number  of  pixels  N  in  a  given 
area.  Although  we  applied  this  technique  primarily  to  EUV / Skylab  spectroheliograms,  it 
is  readily  applicable  to  different  data  sets,  as  long  as  the  difference  in  time  between  two 
consecutive  scans  is  such  that  the  drift  due  to  solar  rotation  within  this  time  interval  is 
much  smaller  than  the  pixel  size  in  the  scans. 
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2.2.3  Image  Processing  Tools 

We  describe  here  two  different  image  processing  tools  that  we  have  found  to  be  powerful 
in  improving  the  spatial  resolution  and  enhancing  the  fine  scale  structures  in  the  images. 
These  are  the  Blind  Iterative  Deconvolution  Algorithm  (BID),  and  the  Image  Enhancement 
Algorithm  (IEA). 

BID  is  a  technique  that  allows  one  to  recover  simultaneously,  from  an  image  with  degraded 
resolution,  the  point  spread  function  (PSF)  and  the  high  angular  resolution  image  through  an 
iterative  process.  It  was  originally  developed  for  the  correction  of  the  effects  of  atmospheric 
turbulence  in  optical  images  (Ayers  and  Dainty  1988)  and  was  modified  for  astronomical  and 
solar  images  (Nisenson  et  al  1989,  Karovska  and  Habbal  1991).  When  applied  to  high  signal- 
to-noise  images,  degraded  by  an  unknown  linear  and  stationary  function,  this  algorithm  finds 
a  pair  of  functions  corresponding  to  the  degrading  function  and  the  undegraded  image,  which 
when  convolved  reproduce  the  input  image  within  given  physical  constraints.  The  first  step 
in  the  iterative  process  is  to  deconvolve  the  degraded  image  with  an  estimate  of  the  PSF  using 
signal-to-noise  ratio  constraints.  The  result  of  the  deconvolution  is  transformed  back  to  the 
image  space,  and  several  physical  constraints,  including  positivity  and  limited  domain  of  the 
image,  are  applied.  Once  the  new  estimate  of  the  PSF  is  obtained,  it  is  used  to  deconvolve 
the  input  image  in  the  next  iteration.  The  iteration  continues  until  some  concluding  criterion 
is  reached,  such  as  image  quality  criterion  or  a  limit  oh  the  number  of  cycles.  Examples  of 
recovered  PSF’s  for  space  observations  from  Skylab  and  Yohkoh  are  shown  in  Figure  lb  and 
lc.  In  both  examples,  the  starting  guess  for  the  PSF  was  a  Gaussian  with  random  noise  and 
half  power  width  radius  of  5  pixels  (Figure  la).  Substantial  sharpening  of  the  structures  in 
the  images  occurred  after  200  iterations,  as  seen  in  the  example  of  the  Yohkoh  image  in  Figure 
Id  and  le,  and  the  iterative  process  was  stopped  when  the  reconstructed  image  remained 
effectively  unchanged  in  further  cycles  of  the  BID.  The  elongation  in  the  reconstructed  PSF 
from  the  Yohkoh  data,  is  probably  due  to  the  off-axis  response  of  the  telescope  mirror.  We 
note  that  there  is  evidence  of  ellipticity  in  the  instrumental  PSF  in  the  measurements  made 
before  the  flight  (see  Maartens  et  al  1994). 

To  explore  the  structure  of  the  solar  atmosphere  at  the  limb,  from  the  chromosphere  to  the 
lower  corona,  we  developed  a  version  of  the  image  enhancement  algorithm  (IEA)  originally 
proposed  by  Koutchmy  et  al  (1988)  (Karovska,  Blundell  and  Habbal  1994).  This  algorithm 
detects  the  location  of  a  significant  change  in  grey  level  in  the  image;  thus,  it  enables  us  to 
obtain  information  regarding  the  spatial  characteristics  of  the  small  scale  structures  despite 
the  presence  of  noise  and  the  high  dynamic  range  between  the  small  and  large  scale  struc¬ 
tures.  This  algorithm  calculates  the  maxima  of  the  multidirectional  second  derivatives  of 
the  intensities  in  a  given  image  using  the  following  nonlinear  operator  (see  Koutchmy  et  al 
1988): 

M(I)  =  aMS[l)  +  (3E2(M8(S2(I)))  +  7E4(M8(52(52(/))))  +  ...  (4) 

The  operator  M8  calculates  the  maximum  of  the  second  derivatives  in  eight  directions,  S2  is 
an  averaging  operator,  E2  and  E4  expand  the  image  with  factors  of  two  and  four,  respectively, 
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Figure  1:  (a)  Starting  PSF  (random  Gaussian),  (b)  reconstructed  PSF  for  a  C  III  Skylab 
image,  (c)  reconstructed  PSF  for  Yohkoh  image  in  (d),  (d)  Yohkoh  image  taken  on  17  April 
1993  of  a  compact  bright  source  at  the  solar  limb,  and  (e)  reconstructed  Yohkoh  image.  The 
lowest  contour  levels  are  at  10%  of  the  maximum,  i.e.  at  the  noise  level. 
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and  I(i,j)  is  the  intensity  in  a  given  pixel  (i,j).  Several  levels  of  block-average  smoothing  may 
be  applied  prior  to  computing  the  derivative,  depending  on  the  signal-to-noise  ratio  in  the 
image,  and  the  characteristic  scale  of  the  structures  in  the  image  that  one  would  like  to 
enhance.  The  coefficients  a,  (3,  and  7  are  weights  chosen  to  select  different  levels  of  block- 
average  smoothing. 


2.3  Results 

2.3.1  Temperature  profile  in  coronal  holes 

Coronal  holes  are  believed  to  be  the  source  region  of  fast  solar  wind  streams  (see,  for  example, 
reviews  in  Zirker  1977).  The  temperature  profile  in  the  inner  corona  is  one  of  the  indicators 
of  the  coronal  heating  process  within  these  regions.  It  also  plays  an  important  role  in 
determining  the  flow  properties  of  the  solar  wind  at  1  AU  (e.g.,  Parker  1958,  Leer  and 
Holzer  1991).  No  single  instrument  to  date  has  measured  directly  the  temperature  profile 
over  a  heliocentric  distance  exceeding  a  fraction  of  a  solar  radius.  Estimates  of  this  profile 
are  derived  from  a  combination  of  inferences,  such  as  shown  in  Figure  2. 


Figure  2:  Temperature  gradient  in  the  inner  corona  as  derived  from  limb  (filled  circles, 
and  open  triangles),  disk  (open  circles  with  filled  centers)  and  in  situ  observations  (triangle 
with  filled  circle).  Arrows  indicate  the  ranges  derived  from  observations.  Hatched  arrows 
represent  measurements  clearly  made  at  and/or  outside  boundaries  of  coronal  holes.  (From 
Habbal,  Esser  and  Arndt  1993). 


Using  Mg  X,  Ne  VII  and  O  VI  line  intensities,  measured  with  the  EUV / Skylab  experiment, 
off  the  solar  limb  in  a  polar  coronal  hole,  Habbal,  Esser  and  Arndt  (1993)  applied  the 
line  ratio  technique  to  derive  the  temperature  below  1.1  Rs.  Although  formed  at  much 
cooler  temperatures,  the  O  VI  line  can  be  used  since  it  has  a  very  broad  response  function  to 
temperatures  reaching  106  K.  The  temperature  across  the  coronal  hole  thus  inferred  is  shown 
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Figure  3:  Temperature  distribution  across  a  south  polar  coronal  hole  and  neighboring  quiet 
regions  as  inferred  from  ratios  of  line  intensities  of  Mg  X,  Ne  VII  and  O  VI.  The  upper  panel 
shows  the  composite  image  in  Mg  X,  where  the  coronal  hole  is  distinguished  by  the  very  low 
emission  with  the  exception  of  bright  points.  The  arrows  indicate  the  radial  directions  along 
which  temperatures  were  inferred.  (From  Habbal,  Esser  and  Arndt  1993). 
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in  Figure  3.  The  top  panel  shows  the  composite  made  from  a  sequence  of  Mg  X  images 
spanning  a  south  polar  coronal  hole  and  its  boundaries.  The  arrows  indicate  the  radial 
directions  along  which  the  line  intensities  were  measured,  and  the  temperature  subsequently 
inferred.  Since  the  intensity  of  the  lines  decreases  very  rapidly  with  heliocentric  distance, 
the  data  are  reliable  only  out  to  1.07  Rs.  The  error  bars  take  into  account  the  uncertainties 
entering  such  an  inference,  as  well  as  the  average  temperature  range  within  the  field  of  view, 
along  a  given  radial  direction. 

A  similar  study  was  carried  out  with  more  recent  ground-based  measurements  of  the  Fe 
X  6374  A  and  Fe  XIV  5303  A  line  intensities  made  at  NSO/SP.  These  lines  have  a  peak 
temperature  of  formation  of  106  K  and  1.8  x  106  K  respectively.  Although  data  were  taken 
for  all  latitudes,  4>,  we  present  here  the  results  for  the  south  polar  coronal  hole  (Figure  4). 
The  boundaries  of  the  coronal  hole  are  very  sharply  defined  by  the  gradient  of  the  Fe  XIV 
intensity  and  were  found  to  match  very  closely  those  determined  from  the  Yohkoh  soft  x-ray 
observations.  Although  made  at  different  heliocentric  distances,  and  different  times,  the 
trend  of  the  temperature  profile  across  the  coronal  hole  is  very  similar  to  that  derived  from 
the  EUV  data  and  shown  in  Figure  3.  It  is  clear,  however,  that  the  use  of  this  combination 
of  iron  lines  fails  in  coronal  holes  when  the  Fe  XIV  emission  is  absent  (see,  for  example,  the 
intensity  of  Fe  XIV  for  150°  <  <fi  <  162°,  at  1.1  Rs,  Figure  4).  When  present,  a  significant 
fraction  of  the  Fe  XIV  emission  is  contributed  by  the  hot  material  veiling  the  coronal  hole 
and  lying  along  the  line  of  sight.  Hence,  the  temperature  inferred  from  these  two  lines  should 
be  regarded  as  an  upper  limit  to  the  coronal  hole  temperature. 

These  two  examples  show  that  temperature  inferences  depend  on  the  choice  of  spectral  lines 
used,  and  that  there  is  a  marked  difference  between  the  temperature  inferred  in  coronal  holes 
and  its  boundaries.  The  larger  temperature  inferred  for  the  boundaries  of  coronal  holes  is  a 
good  indicator  of  the  potential  effect  of  the  line  of  sight  contamination  if  we  envision  that 
such  material  could  also  be  lying  along  our  line  of  sight,  except  in  rare  cases  when  the  extent 
of  the  coronal  hole  along  the  line  of  sight  is  very  large.  This  extent  can  in  fact  be  determined 
from  full  disk  soft  x-ray  observations,  such  as  from  Yohkoh. 

2.3.2  Limits  on  the  helium  abundance  in  coronal  holes 

Another  important  parameter  for  understanding  the  properties  of  the  solar  wind  flow  is  the 
helium  abundance  in  the  inner  corona,  a ,  defined  as  the  ratio  of  helium  to  proton  number 
densities.  A  number  of  theoretical  studies  have  shown  that  the  helium  abundance  in  the  inner 
corona  could  be  significant,  i.e.  »  0.1-0. 2,  and  that  fully  ionized  helium,  or  at  particles, 
could  play  an  important  role  in  the  flow  properties  of  electrons  and  protons  in  the  solar 
wind  (see,  e.g.  Geiss  et  al  1970,  Biirgi  and  Geiss  1986,  Biirgi  1992,  Leer,  Holzer  and  Shoub 
1992).  Observationally,  the  helium  abundance  is  very  hard  to  determine  due  to  the  difficulty 
of  measuring  helium  lines  in  the  corona  (see  discussion  in  Leer,  Holzer  and  Shoub  1992). 
Recently,  we  presented  a  simple  technique  for  inferring  limits  on  the  helium  abundance  in 
the  inner  corona  from  knowledge  of  the  density,  temperature  and  their  respective  gradients 
in  that  region  (Habbal  and  Esser  1994).  We  also  showed  that  in  the  absence  of  density 


Figure  4:  Plot  of  Fe  X  6374  A  intensity,  Fe  XIV  5303  A  intensity,  and  line  ratio  temperature, 
T ,  as  inferred  from  these  two  lines,  versus  latitude  <j),  for  a  south  polar  coronal  hole  observed 
on  10  April  1993.  Measurements  are  shown  at  heliocentric  distances  of  1.1  (solid  line),  1.12 
(short  dashes),  1.14  (dotted  line),  1.16  (long  dashes),  1.18  (dash-dotted  line),  and  1.20  Rs 
(long-short  dashes). 
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measurements,  these  can  be  substituted  by  line  intensity  measurements  and  their  gradients 
in  the  inner  corona.  Using  the  same  EU V / Skylab  data  as  shown  in  Figure  3,  we  inferred 
a  lower  limit  to  a  of  0.3  at  the  base  of  coronal  holes.  For  radial  distances  at  which  the 
temperature  reaches  a  maximum  and  then  starts  to  decrease,  an  upper  limit  to  a  can  be 
derived.  The  example  used  was  for  a  coronal  hole  but  is  readily  applicable  to  different 
coronal  structures.  Thus,  by  applying  this  technique,  inferences  of  limits  on  the  helium 
abundance  in  the  inner  corona  can  be  made  with  data  readily  available  and  to  be  acquired 
in  the  future  either  with  space-  and/or  ground-based  observations. 

2.3.3  Morphology  of  the  magnetic  field  structure  in  coronal  holes 

Knowledge  of  the  plasma  parameters  characteristic  of  coronal  holes  is  but  one  aspect  of  the 
exploration  of  the  coronal  heating  mechanism  (s).  Another  aspect,  often  overlooked,  is  the 
fine  scale  morphology  of  the  magnetic  field  itself.  Macrospicules  and  polar  plumes  are  the 
dominant  structures  in  coronal  holes.  As  seen  in  EU V/ Skylab  observations  made  simulta¬ 
neously  at  different  wavelengths  off  the  limb,  these  structures  are  prominent  at  different 
temperatures  and  have  different  lifetimes  (see,  e.g.  Habbal  1992).  Macrospicules  are  most 
prominent  in  lines  formed  between  104  K  and  105  K.  They  appear  as  columns  of  chromo¬ 
spheric  material,  reminiscent  of  giant  spicules,  extending  anywhere  from  10  to  60”  above  the 
solar  limb.  Their  characteristic  lifetime  ranges  from  3  to  45  min  (Dere  et  al  1989).  Extreme 
ultraviolet  (EUV)  and  radio  observations  suggest  that  they  have  cool  104  K  plasma  cores, 
with  a  thin  outer  sheath  not  exceeding  3  x  105  K  (Withbroe  et  al  1976,  Habbal  and  Gonza¬ 
lez  1991).  The  macrospicules  shown  in  Figure  5,  which  are  the  largest  EUV  macrospicules 
observed  to  date,  have  a  spatial  frequency  of  occurrence  of  approximately  1.5’,  very  simi¬ 
lar  to  that  of  polar  plumes  also  present  in  these  limb  observations  (see  Karovska,  Blundell 
and  Habbal  1994).  While  visible  only  in  lines  formed  above  106  K,  polar  plumes  define  the 
extended  structure  of  coronal  holes  as  seen  in  eclipse  photographs  (e.g.  Koutchmy  1977). 
Recent  work  by  Thieme  et  al  (1989)  suggests  that  remnants  of  polar  plumes  are  observed  in 
interplanetary  space. 

Application  of  the  image  enhancement  algorithm,  IEA,  to  EUV  limb  observations  of  a  north 
polar  coronal  hole  made  with  Skylab,  exposed  the  spicules  at  the  limb  (Figure  5,  bottom 
panel.)  The  substructure  of  the  spicules  and  macrospicules  as  exposed  by  the  application 
of  IEA  is  shown  in  Figure  5.  While  the  limb  appears  as  a  jagged  boundary  in  the  original 
image,  the  processed  image  revealed,  for  the  first  time  in  the  EUV,  a  forest  of  spicule¬ 
like  structures,  with  a  characteristic  spatial  separation  of  10”  to  15”  (Karovska,  Blundell 
and  Habbal  1994).  It  also  revealed  the  existence  of  substructures  within  macrospicules  and 
plumes.  The  processed  images  showed  that  macrospicules  consist  of  column  like  structures, 
with  a  diameter  of  about  5”  -  the  resolution  limit  of  the  telescope,  often  with  a  hook-like 
shaped  tip.  The  processing  also  exposed  low  lying  arches  (about  5”-10”  in  height)  which  were 
commonly  found  to  be  associated  with  their  bases.  The  enhanced  images  of  polar  plumes, 
as  observed  in  O  VI  and  Mg  X,  revealed  the  presence  of  filamentary  substructures  which 
extended  to  at  least  4’,  the  edge  of  the  field  of  view,  with  the  same  characteristic  spatial 
separation  of  10” -15”.  The  enhanced  images  also  showed  that  spicules,  and  the  substructures 
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Figure  5:  Structures  observed  at  the  solar  limb  in  C  III  977  A  in  a  south  polar  region. 
Unprocessed  image  with  contrast  that  emphasizes  macrospicules  (top  panel),  unprocessed 
image  with  contrast  that  emphasizes  the  region  at  the  limb  (middle  panel),  and  image 
processed  with  IEA  (bottom  panel).  (From  Karovska,  Blundell  and  Habbal  1994). 
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within  macrospicules  and  polar  plumes,  observed  at  different  temperatures,  are  not  cospatial. 

2.3.4  Morphology  of  the  magnetic  field  structure  in  active  regions 

The  discrete  nature  of  the  coronal  emission  is  also  a  common  characteristic  of  closed  mag¬ 
netic  structures  such  as  bright  points,  active  and  quiet  regions  (quiet  regions  are  defined  as 
regions  of  diffuse  emission  in  large  scale  closed  magnetic  structures).  Application  of  IEA  to 
simultaneous  multiwavelength  Skylab  observations  showed  that  the  diffuse  emission  in  active 
regions  is  resolved  into  distinct  loops  (Arndt,  Habbal  and  Karovska  1994).  This  technique 
also  confirmed  that  distinct  loops  at  different  temperatures,  such  as  Si  XII  and  Ne  VII 
shown  in  Figure  6,  are  not  cospatial  but  coexist  within  the  same  active  region  complex,  a 
characteristic  also  found  in  the  miniature  loops  forming  coronal  bright  points  (e.g.,  Habbal 
and  Grace  1991).  The  image  enhancement  also  showed  that  the  basic  spatial  scale  in  their 
substructures  is  the  same  in  different  temperature  loops,  a  finding  which  is  not  obvious  in 
the  unprocessed  images.  Similar  results,  i.e.  the  lack  of  cospatiality  and  discrete  nature  of 
the  loops,  were  found  when  IEA  was  applied  to  ground-based  observations  of  Fe  X  and  Fe 
XIV  of  a  quiet  sun  region  observed  off  the  limb  (Figure  7). 


2.3.5  Characteristics  of  the  temporal  variability  of  the  emission 

The  application  of  IEA  to  a  time  sequence  of  the  macrospicule  observations  shown  in  Figure  5, 
showed  the  continuous  interaction  between  their  5”-10”  substructures,  with  the  occasional 
ejection  of  plasma  following  a  ballistic  trajectory  (see  Karovska  and  Habbal  1994).  The 
processed  images  also  revealed  evidence  of  some  pinching  off,  probably  leading  to  the  ejection 
of  plasma  (see,  e.g.  Figure  8  at  23:07:36  UT),  as  had  been  reported  by  Habbal  and  Gonzalez 
(1991)  in  radio  observations  at  4.8  GHz.  The  processed  images  showed  for  the  first  time 
that  the  base  of  macrospicules  is  formed  of  arches  from  which  the  twisted  substructures  that 
form  the  macrospicules  originate  (as  seen  by  comparing  frames  22:56:07  and  22:59:14  UT  in 
Figure  8).  The  morphology  of  the  substructures  within  the  macrospicules  changes  as  well. 
In  fact,  some  of  the  most  dramatic  changes  occur  at  the  base  of  the  macrospicules  where 
the  small  arch-like  structures  reorganize  their  morphology  within  one  minute,  the  shortest 
temporal  resolution  available  in  these  data. 

Applying  the  image  difference  technique  and  image  processing  tools,  to  an  active  region  com¬ 
plex,  we  found  that  the  variable  emission  from  this  complex  is  distinctly  localized  around 
the  footpoints,  as  seen  by  the  white  contours  in  Figure  6  (see  Arndt,  Habbal  and  Karovska 
1994).  In  coronal  holes  and  quiet  regions,  this  variability  is  found  to  be  localized  at  the 
intersection  of  network  cells,  and  in  bright  points  (see  Habbal  and  Grace  1991).  The  spatial 
density  of  the  variability,  6n/N  (as  defined  in  section  2.2.2),  has  a  very  characteristic  tem¬ 
perature  dependence  with  a  peak  at  105K,  in  either  coronal  holes,  quiet  or  active  regions,  as 
shown  in  Figure  9  (Habbal,  Arndt,  and  Karovska  1994).  We  also  point  out  that  the  vertical 
bars,  which  show  the  extent  of  this  measure  as  a  function  of  time,  do  not  extend  by  more 
than  15%.  Another  striking  characteristic  is  that  the  difference  in  the  magnitude  of  Sn/N, 


Si  XII 


Ne  VII 


Figure  6:  Unprocessed  (top  panels)  and  IEA  processed  images  (bottom  panels)  of  an  active 
region  complex  observed  simultaneously  in  Ne  VII  and  Si  XII.The  white  contours  in  the 
lower  panels  indicate  the  location  of  the  significant  variability,  as  determined  by  condition 
(3).  (From  Arndt,  Habbal  and  Karovska  1994.) 
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Figure  7:  Unprocessed  (left  panel)  and  processed  (right  panel)  images  of  a  quiet  region 
observed  in  Fe  XIV  5303  A  at  NSO/SP. 

between  all  three  regions  does  no  exceed  10%.  In  active  regions,  the  magnitude  of  5n/N 
seems  to  be  independent  of  any  flaring  activity,  as  shown  in  Figure  10a.  A  flare  is  evident 
only  in  the  relative  magnitude  change  of  the  intensity  (Figure  10b).  These  results  point 
to  common  characteristic  properties  of  the  underlying  small  scale  magnetic  field  which  are 
independent  of  the  overlying  large  scale  magnetic  structures,  whether  they  are  open  field 
lines  such  as  in  coronal  holes,  or  closed  field  lines  such  as  in  active  or  quiet  regions. 

2.3.6  Coronal  mass  ejections  and  radio  noise  storms 

Acquisition  and  analysis  of  both  disk  and  limb  observations  were  made  to  examine  the 
changing  conditions  of  the  low  solar  corona,  below  1.5  Rs,  preceding  and  following  a  coronal 
mass  ejection  observed  on  the  west  limb  on  12  April  1993.  The  disk  observations  comprised  90 
cm  (333  MHz)  radio  and  daily  Yohkoh  soft  X-ray  measurements,  while  the  limb  observations 
included  measurements  of  emission  from  X-rays,  Fe  X  637.4  nm  and  Fe  XIV  530.3  nm  coronal 
lines,  and  Thomson  scattered  white-light.  The  analysis  of  the  disk  and  limb  observations 
showed  that  throughout  the  three  days  of  consecutive  observations,  the  occurrence  and 
persistence  of  nonthermal  emission  at  90  cm  —  also  known  as  type  I  noise  storm  emission 
—  was  associated  with  the  large-scale  coronal  magnetic  structures  where  the  coronal  mass 
ejection  eventually  occurred  (see  Figure  11).  However,  off-limb  observations  of  the  corona 
gave  no  indication  of  changes,  either  preceding  or  following  the  event,  other  than  a  subsequent 
flare  observed  in  X-rays,  that  were  markedly  different  from  changes  in  other  neighboring 
structures  (Figures  12  and  13).  The  analysis  of  this  novel  combination  of  data  supports  the 
current  view  that  coronal  mass  ejections  are  independent  of  the  “classical”  solar  activity  in 
the  low  corona.  The  radio  observations,  on  the  other  hand,  are  a  strong  advocate  for  the 
connection  between  noise  storms  and  coronal  mass  ejections,  without  implying,  however, 


Figure  8:  A  time  sequence  of  a  section  of  Figure  5,  showing  the  evolution  of  the  central 
macrospicule  as  evidenced  in  the  processed  images.  The  width  of  each  panel  is  75” .  (From 
Karovska  and  Habbal  1994). 
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Figure  9:  ^  versus  temperature  for  an  active  region  (solid  line),  a  quiet  region  (dashed  line) 
and  a  coronal  hole  (dotted  line).  (From  Habbal,  Arndt  and  Karovska  1994). 


Figure  10:  a)  ^  versus  temperature  for  the  regions  defined  by  the  boxes  in  Figure  5.  The 
dashed  curve  is  for  the  flaring  region,  (b)  ~  versus  temperature  for  the  same  sections  of  the 
active  region  as  (a).  The  values  of  ~  for  the  flare  are  given  only  for  Ne  VII  and  Si  XII,  the 
wavelengths  at  which  it  was  recorded,  with  the  vertical  bars  indicating  the  spread  of  this 
measure  as  a  function  of  time.  (From  Habbal,  Arndt  and  Karovska  1994). 
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Figure  12:  Profiles  of  the  Fe  X  637.4  nm  (dash-dot),  Fe  XIV  530.3  nm  (long  dash),  X-ray 
(thick  dash-dot),  at  1.15  Rs,  and  polarization  brightness  intensities  (dot)  at  1.16  Rs,  as  a 
function  of  position  angle  along  the  west  limb  for  7-15  April  1993.  A  few  data  gaps  existed 
in  the  ground-based  observations.  Fe  X  and  Fe  XIV  are  given  in  erg  cm~2  s-1  sr-1;  X-rays 
are  in  arbitrary  units;  and  pB  in  10~7  B/Bs,  where  Bs  is  the  brightness  of  the  solar  disk. 
[From  Habbal  et  al.  1996] 
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that  the  storm  actually  initiates  the  event. 

3.  Conclusion 

The  main  results  from  the  studies  described  in  this  report  lead  us  to  the  following  conclusions: 

Temperature  inferences  in  different  coronal  structures  rely  on  the  spectral  lines  used.  Their 
judicious  choice  is  particularly  crucial  for  reliable  inferences  in  coronal  holes. 

Limits  on  the  helium  abundance  in  the  inner  corona  can  be  inferred  from  knowledge  of  the 
temperature  and  density,  and  their  gradients  in  that  region. 

There  exists  a  characteristic  spatial  separation  of  10”-15”  between  the  substructures  within 
coronal  holes,  that  is  independent  of  the  temperature  of  the  emitting  plasma,  or  the  large 
scale  overlying  magnetic  field. 

Different  temperature  plasmas  coexist  at  coronal  heights  regardless  of  the  overlying  large 
scale  structure,  such  as  active  regions,  quiet  regions  or  coronal  holes,  and  within  these 
structures. 

The  spatial  distribution  of  the  temporal  variability  of  the  coronal  emission  has  a  very  distinct 
temperature  dependence,  with  a  peak  at  105  K,  and  is  independent  of  the  structure  of  the 
overlying  large  scale  magnetic  field. 

In  regard  to  some  of  the  most  spectacular  forms  of  dynamic  magnetic  events,  namely  corona 
mass  ejections,  nonthermal  emission  as  manifested  in  the  form  of  type  I  noise  storms  observed 
at  90  cm,  can  serve  as  predictors  for  the  occurrence  of  these  events. 
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